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The inclusive large-pT production of a single pion, jet or direct photon, and Drell-Yan processes, 
are considered for proton-proton collisions in the kinematical range expected for the fixed-target 
experiment AFTER, proposed at LHC. For all these processes, predictions are given for the trans¬ 
verse single-spin asymmetry, Av, computed according to a Generalised Parton Model previously 
discussed in the literature and based on TMD factorisation. Comparisons with the results of a 
collinear twist-3 approach, recently presented, are made and discussed. 

PACS numbers: 13.88.-l-e, 13.85.Ni, 13.85.Qk 


I. INTRODUCTION AND FORMALISM 

Transverse Single-Spin Asymmetries (TSSAs), have been abundantly observed in several inclusive proton-proton 
experiments since a long time; when reaching large enough energies and pT values, their understanding from basic 
quark-gluon QCD interactions is a difficult and fascinating task, which has always been one of the major challenges 
for QCD. 

In fact, large TSSAs cannot be generated by the hard elementary processes, because of helicity conservation (in 
the massless limit) typical of QED and QCD interactions; indeed, such asymmetries were expected to vanish at high 
energies. Their persisting must be related to non perturbative properties of the nucleon structure, such as parton 
intrinsic and orbital motion. A true understanding of the origin of TSSAs would allow a deeper understanding of the 
nucleon structure. 

Since the 1990s two different, although somewhat related, approaches have attempted to tackle the problem. 
One is based on the collinear QCD factorisation scheme and involves as basic quantities, which can generate single 
spin dependences, higher- twist quark-gluon-quark correlations in the nucleon as well as higher-twist fragmentation 
correlators. The second approach is based on a physical, although unproven, generalisation of the parton model, with 
the inclusion, in the factorisation scheme, of transverse momentum dependent partonic distribution and fragmentation 
functions (TMDs), which also can gen erate single spin dependences. The twist-3 correlations are related to moments 
of some TMDs. We refer to Refs. [ll-Q: and references therein, for a more detailed account of the two approaches, 
and possible variations, with all relevant citations. Following Ref. 0, we denote by CT-3 the first approach while 
the second one is, as usual, denoted by GPM. 

In this paper we consider TSSAs at the proposed AFTER@LHC experiment, in which high-energy protons extracted 
from the LHC beam would collide on a (polarised) fixed target of protons, with high luminosity. For a description 
of the physics potentiality of this experiment see Ref. [ll| and for the latest technical details and importance for 
TMD studies see, for example. Ref. [l^. Due to its features the AFTER@LHC is an ideal experiment to study 
and understand the origin of SSAs and, in general, the role of QCD interactions in high-energy hadronic collisions; 
AFTER@LHC would be a polarised fixed target experiment with unprecedented high luminosity. 

We recall our formalism by considering the Transverse Single-Spin Asymmetry Ajv, measured in—>■ hX inclusive 
reactions and defined as: 


_ da^ - da^ 
dcrl -b dcrl 


with 


da^’^ 


d^Ph 


( 1 ) 


where t,i are opposite spin orientations perpendicular to the x-z scattering plane, in the pp^ c.m. frame. We define 
the t direction as the -by-axis and the unpolarised proton is moving along the -bi-direction. In such a process the 
only large scale is the transverse momentum px = \ {Ph)x\ of the final hadron. 

In the CPM Am originates mainly from two spin and transverse momentum effects, one introduced by Sivers in 
the partonic distributions [I10, and one by Collins in the parton fragmentation process [T^, being all the other 
effects strongly suppressed by azimuthal phase integrations [l6|. According to the Sivers effect the number density 
of unpolarised quarks q (or gluons) with intrinsic transverse momentum k± inside a transversely polarised proton p^, 
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with three-momentum P and spin polarisation vector 5, can be written as 

fq/pt{x,k±) = fq/p{x,k±) -h i A^fg/p-t{x,k±) S- {Pxk±), (2) 

where x is the proton light-cone momentum fraction carried by the quark, fq/p{x,k±) is the unpolarised TMD 
{k± = |fe_L|) and A^fq/pt{x,k±) is the Sivers function. P = P/\P\ and k± = k±/k± are unit vectors. Notice that 
the Sivers function is most often denoted as fiT{x,k±) [l3|; this notation is related to ours by [l8j | 

^^fq/prix, k±) = f^r^ix, kj_) , (3) 

TJlp 

where rup is the proton mass. 

Similarly, according to the Collins effect the number density of unpolarised hadrons h with transverse momentum 
resulting in the fragmentation of a transversely polarised quark , with three-momentum q and spin polarisation 
vector Sq^ can be written as 

Dqyh{z,P±) = Dh/q{z,p±) -f i A^Dgt/f,{z,p±) Sq- {qx p_l) , (4) 


where z is the parton light-cone momentum fraction carried by the hadron, D^/qiz^p^) is the unpolarised TMD 
(p± = Ip_lI) and A^Dqt/h{z,p±) is the Collins function, q = q/\q\ and = p±/p± are unit vectors. Notice that 
the Collins function is most often denoted as H^'^{z,p±) [13; this notation is related to ours by (l8j | 

H^^z,p±) , (5) 

where is the hadron mass. 

According to the GPM formalism A]\[ can then be written as: 


An = 


[da'^ - do-*'] Sivers + [da^ - do-'^-jcollins 


do'l' -I- do'-l 

The Collins and Sivers contributions were recently studied, respectively in Refs. [3 and [3, and are given by: 


( 6 ) 


[dc'^ do"^] Sivers — ^ ^ 

a,6,c,(i 


dxa dxb dz 2 , .2 


16 TT^ Xa Xb z'^S 


d k±a d^k±b d Px S{pj_ ■ p^) J{pi_) S{s + t + u) 


1 r, 


X A fa/pf{Xa,kxa) COs{(j)a) fb/pixb,kxb) + \M^\ + \M^ 


f0|2 


- ab—¥cd 


Dh/ciz,P±) , (7) 


and 


[do-7 - dcr-*-] Collins = / 


dxa dxb dz 
16 TT^ Xa Xb z'^S 


d^kxa d^k±b d^Pj_ Hp± ' Pc) JiP±) + i + u) 


qa,b,qc,d *' 

X ArqaiXa^k^a) COs{4>a + Pi - <P2 + 4'^) 
X fb/p(^Xb^kj_b) 


( 8 ) 




A^D 


qab^qcd 


h/ql^^^P^)- 


For details and a full explanation of the notations in the above equations we refer to Ref. (TB] (where p^ is denoted 
as k±c)- It suffices to notice here that J{p±) is a kinematical factor, which at 0{p±/Eb) equals 1. The phase factor 
cos(^a) in Eq. (O originates directly from the k± dependence of the Sivers distribution [S ■ {P x k±), Eq. ([2])]. The 
(suppressing) phase factor cos{(j)a + <pi — P 2 + 4^) in Eq. ([8|) originates from the k± dependence of the unintegrated 
transversity distribution Axq, the polarized elementary interaction and the spin-pj^ correlation in the Collins function. 
The explicit expressions of ipi^Lp 2 and (j)^ in terms of the integration variables can be found via Eqs. (60)-(63) in [l^ 
and Eqs. (35)-(42) in (T^j- 

The MO’s are the three independent hard scattering helicity amplitudes describing the lowest order QCD interac¬ 
tions. The sum of their moduli squared is related to the elementary unpolarised cross section that is 


di 


1 1 
167rs2 2 


2=1 


( 9 ) 
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The explicit expressions of the combinations of M?’s which give the QCD dynamics in Eqs. © and dH), can be 
found, for all possible elementary interactions, in Ref. [l^ (see also Ref. [l| for a correction to one of the product of 
amplitudes). The QCD scale is chosen as Q = pt- 

The denominator of Eq. © or (ini) is twice the unpolarised cross section and is given in our TMD factorisation by 
the same expression as in Eq. ©, where one simply replaces the factor cos{(f>a) with 


II. An for single PION, JET AND DIRECT PHOTON PRODUCTION 

We present here our results for A^r, Eq. ([T|), based on our GPM scheme, Eqs. ([H]), © and dH). The TMDs which 
enter in these equations are those extracted from the analysis of Semi Inclusive Deep Inelastic (SIDIS) and e+e" 
data adopting simple factorised forms, which we recall here. For the unpolarised TMD partonic distributions 

and fragmentation functions we have, respectively: 


fq/p{x, k±) = fq/p{x) 


(fc2) =0.25GeV2 


( 10 ) 


and 


e-pi/(pi> 

Dh/q{z,p±) = Di^fq{z) ^^^2 ^ 


The Sivers function is parameterised as 


=o.20GeV^ 




where 


fqlp^ (.X, k±) 2 (x) fq/p{x) h{k±) ^ 

MHx) = <x“^(l - 


) ’ 


with |A^| < 1, and 


cTq’^fq^ 




( 11 ) 


( 12 ) 


(13) 


(14) 


Similarly, the quark transversity distribution, ^T<l{x,k\_\ and the Collins fragmentation function, Df^/qt{z,p±), 
have been parametrized as follows: 


ATq{x,k±) = -Afq{x) [fq/p{x) + Aq{x)] 


o-kl/{kl) 

7r(fc2) 


(15) 


^ Dh/q-f{.Z, PJ_) = 2Af^(z) D/,/q(z) h(p±) ^ , 


(16) 


where Aq{x) is the usual collinear quark helicity distribution. 




AG (x) = Ajx^Gl-a;) " 


dn 

dq Qjq 


(17) 


A7f (z) = Af ^ 


Iq 


(18) 


with < 1, and 


h{p±) = \f2e ^ e 


( 19 ) 
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FIG. 1: Our theoretical estimates for ^jv vs. pr at ^/s = 115 GeV, xp = —0.2 (upper plots) and xf = —0.4 (lower plots) for 
inclusive TT^ and TT^ production in pp^ —>■ ttX processes, computed according to Eqs (El)-® of the text. The contributions 
from the Sivers and the Collins effects are added together. The computation is performed adopting the Sivers and Collins 
functions of Refs. [^ . (SIDIS 1 - KRE, left panels), and of Refs. [^, (SIDIS 2 - DSS, right panels). The overall 
statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following 
the procedure described in Appendix A of Ref. . 


All details concerning the motivations for such a choice, the values of the parameters and their derivation can be 
found in Refs. [20l - l^ . We do not repeat them here, but in the caption of each hgure we will give the corresponding 
references which allow to fix all necessary values. 

We present our results on A at for the process —>• ttX at the expected AFTER@LHC energy (-/s =115 GeV) in 

Figs. [T1I51 Following Refs. 0,1, our results are given for two possible choices of the SIDIS TMDs, and are shown as 
function of pp at two hxed xp values (Fig. [T]), as function of xp at two fixed rapidity y values (Fig. [2|) and as function 
of rapidity at one fixed pp value (Fig.®), xp is the usual Feynman variable defined as 0 :^ 7 ’ = 2pply/s where pi = {Ph)z 
is the z-component of the final hadron momentum. Notice that, in our chosen reference frame, a forward production, 
with respect to the polarised proton, means negative values of xp. The uncertainty bands reflects the uncertainty in 
the determinations of the TMDs and are computed according to the procedure explained in the Appendix of Ref. 211. 
More information can be found in the figure captions. 

Notice that, for both our choices of the Sivers functions, the gluon Sivers distributions are taken to be vanishing, as 
suggested by data [2l|, [ 2 ^. Gluon channels contribute instead to the unpolarised cross sections, in the denominator 
of Eq. ([T]) or (®. For the unpolarised partonic distributions we adopt the GRV98LO PDF set [ 2 ^ and for the 
fragmentation functions the DSS set from Ref. [ 2 ^ and the Kretzer (KRE) set from Ref. ( 2 ^ . 

The analogous results for the single direct photon are shown in Figs. HHHl (where xp = ^/^), and those for 

the single jet production in Figs. [7][n] {xp = 2{p^)z/y/s). In these cases, obviously, there is no fragmentation process 
and only the Sivers effect contributes to Api, with D}^/^{z,pp) simply replaced by 5{z — l)(5^(pj.) in Eq. (® (see 
Ref. 0 for further details). In our leading order treatment the jet coincides with a single final parton. Notice that for 
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FIG. 2: Our theoretical estimates for An vs. xf at ^/s = 115 GeV, y = —1.5 (upper plots) and y = —3.0 (lower plots) for 
inclusive and 7r° production in pp^ —>■ ttX processes, computed according to Eqs (El)-® of the text. The contributions 
from the Sivers and the Collins effects are added together. The computation is performed adopting the Sivers and Collins 
functions of Refs. [^ . (SIDIS 1 - KRE, left panels), and of Refs. [^, (SIDIS 2 - DSS, right panels). The overall 
statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following 
the procedure described in Appendix A of Ref. f2ll| . 


a jet production we have all the same QCD subprocesses which contribute to hadron production, while for a direct 
photon production the basic partonic subprocesses are the Compton scattering g q{q) —> 7 q(q) and the annihilation 
process qq ^ "f g [^ . 


III. An for DRELL-YAN PROCESSES 

Drell-Yan (D-Y) processes are expected to play a crucial role in our understanding of the origin, at the partonic 
level, of TSSAs. For such processes, like for SIDIS processes and contrary to single hadron production, the TMD 
factorisation has been proven to hold, so that there is a general consensus that the Sivers effect should be visible via 
TSSAs in D-Y Not only: the widely accepted interpretation of the QCD origin of TSSAs as final or initial 

state interactions of the scattering partons |33l| leads to the conclusion that the Sivers function has opposite signs in 
SIDIS and D-Y processes [s^. Which remains to be seen. 

Predictions for Sivers Ajy in D-Y and at different possible experiments were given in Ref. [35l |. which we follow here. 

In Ref. ( 3 ^ predictions were given for the p^p —>■ t'^l~X D-Y process in the p^ — p c.m. frame, in which one observes 
the four-momentum q of the hnal t^t~ pair. Notice that q^ = is the large scale in the process, while qp = IqtI is 
the small one. In order to collect data at all azimuthal angles, one defines the weighted spin asymmetry: 

^ lo^ d(j)^ - da^\ sin(^^ - ^g) 

^ 1 lo^ ■*" 


( 20 ) 
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FIG. 3: Our theoretical estimates for An vs. y at y/s = 115 GeV and px = 3 GeV, for inclusive and tt® production in 

pp^ TT X processes, computed according to Eqs dS])-® of the text. The contributions from the Sivers and the Collins effects 
are added together. The com put ation is performed adopting the Sivers and Gollins functions of Refs. [13, [l^l (SIDIS 1 - KRE, 
left panel), and of Refs. [U, (SIDIS 2 - DSS, right panel). The overall statistical uncertainty band, also shown, is the 
envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix A of 
Ref. [^. 


J d(j)-y 

X), / d'^kpi d?kF2 + fcj_2 - Qt) ^^fq/prixi,k±i) fg/pix2,k±2) 

sin 

(pj (ps) 

f d4>-y 

J 2 q^q f d‘^kj_l d?kF2 d‘^{k±l + A;f2 - Qt) fq/p{xi,kFl) fq/p{x2,kj_2) 



( 21 ) 


where 4>-y and (j)s are respectively the azimuthal angle of the pair and of the proton transverse spin and we have 
defined (see Eq. @): 


= A^fg/pr{x,k±) S ■ {P x k±) = fq/p-t{x,k±) - fg/pi{x,ki_) . (22) 

Adopting for the unpolarised TMD and the Sivers function the same expressions as in Eqs. (I10|) and (ll2D - (n3 
allows, at 0{k±/M), an analytical integration of the numerator and denominator of Eq. (I21L resulting in a simple 
expression for the asymmetry jssj. 

Notice that we consider here the pp^ —>■ i^£~X D-Y process in the p — p^ c.m. frame. For such a process the TSSA 
is given by [s^ 


^ XF, M, qr) = ^ 7 *A; -xp, M, qt) . (23) 

Our results for the Sivers asymmetry at AFTER@LHC, obtained following Ref. [s^, Eq. (l23l) and using 

the SIDIS extracted Sivers function reversed in sign, are shown in Fig. 1101 Further details can be found in the captions 
of these figures. 


IV. COMMENTS AND CONCLUSIONS 


Some final comments and further details might help in understanding the importance of the measurements of the 
TSSAs at AFTER@LHC. 


• Most predictions given show clear asymmetries, sufficiently large as to be easily measurable, given the expected 
performance of AFTER@LHC ll| . The uncertainty bands reflects the uncertainty in the extraction of the Sivers 
and transversity functions from SIDIS data, which are focused on small and intermediate x values {x < 0.3); in 
fact the bands grow larger at larger values of |a:F|- 


• The values of A^r found for pion production can be as large as 10% for while they are smaller for 7r°. They 
result from the sum of the Sivers and the Collins effects. The relative importance of the two contributions varies 
according to the kinematical regions and the set of distributions and fragmentation functions adopted. As a 
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FIG. 4: Our theoretical estimates for An vs. pt at ^/s = 115 GeV, xp = —0.2 (upper plots) and xf = —0.4 (lower plots) for 
inclusive photon production in pp^ —>■ 7 X processes, computed according to Eqs Band 0 of the text. Only the Sivers effect 
contributes. The computation is performed adopting the Sivers functions of Ref. |20ll (SIDIS 1, left panels) and of Ref. (2l| 
(SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in 
Appendix A of Ref. [Ij. 


tendency, the contribution from the Sivers effect is larger than the Collins contribution with the SIDIS I - KRE 
set, while the opposite is true for the SIDIS 2 - DSS set. 

The values found here are in agreement, both in sign and qualitative magnitude, with the values found in 
Ref. [13 within the collinear twist-3 (CT-3) approach. 

• The results for single photon production are interesting; they isolate the Sivers effect and our predictions show 

that they can reach values of about 5%, with a reduced uncertainty band. We find positive values of Ayy as the 
relative wei ght of the quark charges leads to a dominance of the u quark and the Sivers functions is 

positive [^.1^. 

Our results, obtained within the GPM, have a similar magnitude as those obtained in Refs. M and Q, within 
the CT-3 approach, but have an opposite sign. Thus, a measurement oi An for a single photon production, 
although difficult, would clearly discriminate between the two approaches. 

• The values oi An for single jet production, which might be interesting as they also have no contribution from 
the Collins effect, turn out to be very small and compatible with zero, due to a strong cancellation between the 
u and d quark contributions. The same result is found in Ref. [I 3 . 

• A measurement of in D-Y processes at AFTER@LHC is a most interesting one. In such a case the 

TMD factorisation has been shown to be valid and the Sivers asymmetry should show the expected sign change 
with respect to SIDIS processes [H, [s^. Our computations. Fig. [TOl predict a clear asymmetry, which can be 
as sizeable as 10%, with a definite sign, even within the uncertainty band. 
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FIG. 5: Our theoretical estimates for An vs. xf at ^/s = 115 GeV, y — —1.5 (upper plots) and y = —3.0 (lower plots) for 
inclusive photon production in pp^ —>■ 7 X processes, computed according to Eqs (151) an d 0 of the text. Only the Sivers effect 
contributes. The computation is performed adopting the Sivers functions of Ref. |20ll (SIDIS 1, left panels) and of Ref. (2l| 
(SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in 
Appendix A of Ref. [Ij. 


Both the results of Ref. [13 and the results of this paper, obtain solid non negligible values for the TSSA An 
measurable at the AFTER@LHC experiment. The two sets of results are based on different approaches, respectively 
the CT-3 and the GPM factorisation schemes. While the magnitude of Ajv is very similar in the two cases, the signs 
can be different; in particular, the TSSA for a direct photon production, pp^ —> 7 A, has opposite signs in the two 
schemes. 

In this paper we have also considered azimuthal asymmetries in polarised D-Y processes, related to the Sivers 
effect. As explained above, in this case, due to the presence of a large and a small scale, like in SIDIS, the TMD 
factorisation is valid, with the expectation of an opposite sign of the Sivers function in SIDIS and D-Y processes. Also 
this prediction can be checked at AFTER@LHC. 
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FIG. 6: Our theoretical estimates for ^jv vs. y at y/s = 115 GeV and pr ~ 3 GeV, for inclusive photon production in 
—>■ y X processes, computed according to Eo s (l6l) -dH) of the text. Only the Sivers effects contributes. The computation is 
performed adopting the Sivers functions of Ref. [2^ (SIDIS 1, left panel) and of Ref. [2^ (SIDIS 2, right panel). The overall 
statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following 
the procedure described in Appendix A of Ref. . 
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FIG. 7: Our theoretical estimates for Ajv vs. xf at y/s = 115 GeV, xf = —0.2 (upper plots) and xf = —0.4 (lower plots) 
for inclusive single jet production in —>■ jet A processes, computed according to Eqs and © of the text. Only the 

Sivers effect contributes. The computation is performed adopting the Sivers functions of Ref. [2^ (SIDIS 1, left panels) and 
of Ref. (SIDIS 2, right panels}. The overall statistical uncertainty band, also shown, is obtained following the procedure 
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FIG. 8: Our theoretical estimates for An vs. xf at ^/s = 115 GeV, y — —1.5 (upper plots) and y = —3.0 (lower plots) for 
inclusive single jet production in pp^ —>■ jet X processes, computed according to Eqs ® and 0 of the text. Only the Sivers 
effect contributes. The computation is performed adopting the Sivers functions of Ref. (SIDIS 1, left panels) and of Ref. |2l| 
(SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in 
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FIG. 10: Our theoretical estimates for in D-Y processes as expected at AFTER@LHC. Our results are presented as 

function of M (upper plots), xf (middle plots) and x of the quark inside the polarised proton, (lower plots). The other 
kinematical variables are either fixed or integrated, as indicated in each figure. They are computed according to Ref and 
Eq. (I23II . adopting the Sivers functions of Ref. (SIDIS 1, left panels) and of Ref. (SIDIS 2, right panels), reversed 
in si gn. The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix A of 
Ref. [in. 
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